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O \ Electron tunneling experiments are used to probe Coulomb correlation effects 

in the single-particle density-of-states (DOS) of boron-doped silicon crystals 

'^ \ near the critical density Uc of the metal- insulator transition (MIT). At low 

energies {e < 0.5 meV), a DOS measurement distinguishes between insulating 
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^D ' and metallic samples with densities 10 to 15 % on either side of Uc- However, 
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^^ ' at higher energies (~ 1 meV < e < 50 meV) the DOS of both insulators 
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Cj I and metals show a common behavior, increasing with energy as e™" where m 

'^ ' is roughly 0.5. The observed characteristics of the DOS can be understood 
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using a classical treatment of Coulomb interactions combined with a phe- 
nomenological scaling ansatz to describe the length-scale dependence of the 
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$H ' dielectric constant as the MIT is approached from the insulating side. 



71.23.-k, 71.30.-hh, 71.45.Gm 



Typeset using REVT^ 



I. INTRODUCTION 

Since 1979 the best available description of the disorder-driven metal-insulator transition 
(MIT) has been based on the non-interacting scaling theory of Abrahams, et al. |]l| Although 
it treats thoroughly the localizing effects of electron scattering off static disorder, this theory 
is incomplete because it neglects electron-electron interactions. In particular, below the 
critical density ric of the MIT the vanishing carrier mobility means that Coulomb correlations 
are strong enough to warrant treatment on equal footing with the disorder. The importance 
of interactions is illustrated by several well-known, apparently anomalous phenomena, |0] 
such as the recent discovery by Kravchenko, et al. [^,^ of a "forbidden" metallic state in a 
two-dimensional (2D) electron system. In a disordered metal, electron-electron interactions 
lead to a singular negative correction to the single-particle density-of-states (DOS) near the 
Fermi level, f^. Deep into the insulating side, it is well established that Coulomb correlations 
cause a Coulomb gap in the DOS near the Fermi level, which changes the temperature 
dependence of the DC electrical conductivity at very low temperatures. PJ^ How these two 
renormalizations of the DOS evolve in the critical region of the MIT and match each other 
at He is one of the most challenging and long-standing questions of solid-state physics. 

Combining both disorder and Coulomb interactions into a unified scaling description of 
disordered metals near the MIT was attempted by McMillan [^ and by Gefen and Imry [Q . 



These models were criticized ||10l for the use of the single-particle, rather than thermody- 
namic, DOS in describing the charge screening, as discussed in detail in Sect. Ill, and 
have not been widely accepted despite garnering some experimental support. |jll| Later 



renormalization group theories by Finkelstein, [|1^ Castellani, et al, ||13[ and Kirkpatrick 
and Belitz 0] used a Fermi-liquid based approach to describe the diffusion of interacting 
quasiparticles in a disordered metal. These works employed a smooth thermodynamic DOS 
(or compressibility) and started from perturbative treatments of weak disorder. All agree 
with the basic findings of Altshuler and Aronov |Q (AA) that there should be a square-root 
(in 3D) or logarithmic (in 2D) depletion in the single-particle DOS at the Fermi level of a 



disordered metal. However, they failed to generate a well-defined, continuous charge local- 
ization transition at nonzero disorder and finite interaction strength, and cannot describe 
the emergence of a Coulomb glass state and the Coulomb gap as a system crosses over into 
the insulating state. 

Tunneling experiments have observed AA-like depletions in the DOS of a variety of dis- 
ordered metals, including amorphous metal- semiconductor alloys, ||l5| , p!6| , |TT| , doped semi- 



conductors, [13] and granular metals. |I^ By contrast, with the exception of the singular 
sodium tungsten bronze system, [^ the existence of the Coulomb gap in localized insulators 
has until recently been only indirectly inferred from activation fits to DC conductivity pOl or 



relaxation measurements. |21l| Only in the last few years have quantitative tunneling spectro- 



scopic observations of the Coulomb gap been made in 3D p2| , p3| and in 2D p4| non-metallic 
semiconductors. 

Most previous work on interaction effects has emphasized the metallic side; much less 
corresponding effort, experimental or theoretical, exists to describe the insulating side near 
the critical region of the MIT. In this paper we concentrate on studying the Coulomb gap 
in the DOS that occurs in 3D disordered insulators close to the MIT. Results of tunneling 
measurements of the DOS in boron-doped silicon are presented over a much larger range of 
energies than in Ref p2| , |23|| . The tunneling DOS spectra show that metals and insulators 
can be distinguished by the low-energy characteristics of the DOS {i.e. square-root cusp vs. 
parabolic Coulomb gap), but that metals and insulators share a common, roughly square- 
root high-energy DOS behavior that is approximately independent of dopant densities close 
to ric. Of course, a truly complete microscopic theory of the MIT with interactions should 
cover the continuous crossover from insulating to metallic state. We present a much more 
modest and simple scaling approach to the MIT, which emphasizes the insulating side and 
is based on an extension of the Efros-Shklovskii [^] arguments. Using the idea of a dielectric 
constant with spatial dispersion, [§,|| we argue that at the point of the MIT a power law 
depletion of the DOS near the Fermi level is still the ES Coulomb gap, only modified by this 
dispersion. The scaling ansatz describes well the main features of the experimental data. 



This paper is organized as follows: Section II presents experimental details and data, 
Section III develops the scaling description, and Section IV analyzes of the data with respect 
to the scaling model. 

II. TUNNELING MEASUREMENTS OF THE COULOMB GAP 

When a conductor is separated from a conventional metal by a rectangular potential 
barrier high enough to prevent classical current flow but thin enough to permit quantum 
tunneling, the tunneling conductance at temperature T, G{V,T) = dl/dV where I is the 
tunneling current and V is the voltage bias across the junction, is given by: p5 

G{V) fN(E) af(e-eV,T) 



Go J iV„ aev <" <") 

where Go is the conductance in the non- interacting case, N{e, T)/Nq is the interacting single- 
particle DOS relative to the non-interacting value A'o, and / is the Fermi function. We take 
the zero of energy at ep = 0. The highest bias used in this experiment is 50 meV, which is 
much smaller than the several eV height of the Si02 barrier, so that the barrier transmission 
coefficient is taken to be independent of bias. In nearly all tunneling measurements, the 
normalizing conductance Gq is taken at a relatively high voltage bias, where G{y) is either 
constant or only slowly varying on the energy scale of the spectral features of interest. 
The normalized conductance then gives the ratio N{eV,T)/NQ, thermally broadened by 
convolution with —df/d{eV). At sufficiently low temperature, when ksT is much smaller 
than the energy scale of characteristic variations in N{e), ordinary thermal broadening can 
be neglected so that G{V,T)/Go is directly proportional to A^(eV,T)/A^o- Note that, unlike 
the non-interacting case, N{eV, T) can have a non-trivial intrinsic temperature dependence 
separate from ordinary thermal broadening. 

The Si:B crystals used were characterized by measuring their room-temperature resistivi- 
ties p and resistivity ratios (RRs) p(4.2 K)/p(300 K). Dopant densities n were obtained using 
the calibration of Thurber, et al. |2^ and the data of Dai, et al. [^ to translate between 



measured RRs and dopant densities. Samples used had RRs of 2.3 to 18, corresponding 
to n/rif. of 110 % to 81 %, respectively, where we take |2^ ra^ = 4.0 x 10^*^ cm~^. Details 



of the DC variable-range hopping (VRH) conductivity in many of these samples have been 
published previously. p^j28[] The static dielectric constants of several insulating crystals 
were measured using a capacitance method similar to that described by Castner. IQ These 
results are shown in Fig. |l[ To compare values with the literature, on an 86 % sample we 
obtained k/kq = 8.5 ± 1, where kq = 11-7 is the dielectric constant of the host Si lattice. 
This value is in reasonable agreement with published values |^ on Si:P and Si:As at similar 
values of n/ric- The curve in Fig. |l] is a fit of the data to K{n) ~ (1 — n/ric)^'^, which yields 
C~0.71. 

Tunnel junctions were formed by creating a very thin (estimated 15 to 30 A) layer 
of Si02 as a tunnel barrier on the Si:B crystals, and using an Al film as the counter- 
electrode. Details of the fabrication of these metal-oxide-semiconductor structures were 
described previously. The major modification from the prior description was the use of an 
ultraviolet ozone method in place of a chemical oxide growth to clean the silicon surfaces and 
produce more reproducible ultra-thin Si02 tunnel barriers. Si02 thicknesses were estimated 
from published calibrations of UV exposure time, temperature, and O2 flow rate. M^ Good 



junctions were those which showed a signature of the Al superconducting energy gap below 
1 K. These devices had junction resistances between 65 to 850 kQ near 1 K. 

Tunneling current- volt age (/ — V) and conductance- voltage {G{V)) traces were taken 
using standard analog methods. Where required, a small (1 kG) magnetic fleld was used 
to suppress the superconductivity in the Al electrode. We shall refer to such data as "zero 
Teslas." Data down to 1.2 K were obtained by suspending the samples stress- free from their 
leads and immersing in a pumped liquid ^He bath. To reach lower temperatures, the samples 
were immersed in the ^He/^He mixture of a dilution refrigerator. Cooling power at 0.1 K 
was measured to be 120 /iW, while no more than 0.1 fiW of power was used to take the 
data. Temperature stability of better than 1 niK at 1.2 K and 0.1 mK at 0.1 K could be 
obtained, which was important to prevent thermal fluctuations from coupling to the highly 



temperature-dependent resistance in the more insulating samples. For samples with n/ric > 
90 %, reliable data could be taken below 0.1 K. For insulating samples with n/ric < 90 
%, the resistance of the Si:B crystal itself rose to exceed 10 % of the resistance across the 
tunnel junction below 0.1 K. Because the crystal acted as a voltage drop in series with the 
junction, quantitative tunneling conductance data below ~ 0.1 K is not reliable in these most 
insulating samples. Therefore, we present only T > 0.5 K data for these samples, where the 
voltage drop across the Si:B crystal is no more than a 1 % correction to the tunnel junction 
conductance. 

Figure § shows low-energy tunneling conductance spectra on six samples with densities 
n/ric = 81 %, 86 %, 90 %, 95 %, 105 %, and 110 %, all in zero Teslas. All data are normalized 
to the conductance value of each sample at +2 mV. The spectrum for the 81 % sample was 
taken at 1.2 K, the spectrum for the 86 % sample was taken at 0.5 K, and the data for the 
other four samples were taken at 0.1 K. Ordinary thermal broadening at the measurement 
temperatures is minor compared to the energy scale of the conductance features, and so has 
not been deconvoluted from the data. However, the interacting DOS of each sample has 
an intrinsic finite temperature dependence that may be more obvious in the data for the 
81 % and 86 % samples because of the higher temperature. In the conductance spectra, 
there is a clear difference between the insulating and metallic samples at low biases, |V^| < 
1 mV. In all the insulating samples, a conductance dip across zero voltage (the Fermi level 
in tunneling measurements) is the signature of the Coulomb gap. For the 81 % and 86 % 
samples, G{V = 0) is small but greater than zero because the measurement temperature 
for these two samples is sufficiently high that there is some intrinsic filling in of the T = 
gap. In all insulating samples, the low-energy conductance spectra fit a power-law form 
N{e) oc \e\P with p = 2.4, 2.2, 2.2, and 2.0 for the 81 %, 86 %, 90 %, and 95 % samples 
respectively, at the (different) measurement temperatures. As n — > ra^ from below, the gap 
sides steepen and the width of the Coulomb gap around zero bias narrows. Although the gap 
"shoulders" are soft and therefore not precisely definable, the approximate full width delta 
of the parabolic-like {p ~ 2) gap is about 2.2 mV, 1.5 mV, 1.1 mV, and 0.6 mV and for the 81 
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%, 86 %, 90 % and 95 % samples respectively. By contrast, both the metallic samples show 



a significantly large nonzero conductance at V = 0, with a sharp dip |^ around zero bias 
that is well fit to a square-root form N{e) = N (0)[1 + (\£\ / 5Y^'^] at low bias. Fits to this form 
give 6 ~ 0.25 meV and 0.4 meV for the 110 % and 105 % samples, respectively. Normalizing 
the metallic data to G(+50 mV), the ratio of zero-bias conductances is G'iio(0)/G'io5(0) = 
1.7, so that the more metallic sample has a larger DOS at the Fermi level, as expected. 

Figure || compares tunneling conductance spectra over an extended bias range (0.03 to 
50 mV) for two metallic and two insulating samples 5 % and 10 % above and below ric, all at 
0.1 K. This log-log-plot clearly reveals a characteristic energy scale of approximately 0.5 to 
1 mV. Below ~ 0.5 mV, an obvious distinction can be made between metallic and insulating 
samples. The low-bias data for the metallic samples approach V^ = with a thermally 
rounded square-root shape and have a nonzero G(y = 0). By contrast, the insulating 
samples show a quasi-parabolic Coulomb gap depletion of the low-energy DOS. The most 
important feature of Fig. |^ is that from ~ 1 mV to 50 mV, the conductance spectra for both 
metals and insulators are essentially indistinguishable. The high-bias tunneling conductance 
common for both insulating and metallic samples follows a functional form G{V) oc V"^ 
where m is between 0.43 to 0.47, with no correlation of m with dopant concentration in the 
range of samples studied. Thus above a characteristic energy a DOS measurement cannot 
differentiate between metallic and insulating states. A similar dependence was reported by 



Hertel, et al. |TT| in tunneling studies of metallic NbSi alloys near the MIT. They found 
that the tunneling conductance showed a square-root cusp at low energies that turned over 
to a slower than square root dependence at higher energies. In NbSi m was measured to be 
closer to 1/3. 

While all the junctions used Si02 barriers of roughly equal thicknesses, the common 
high-bias conductance behavior for metals and insulators is very unlikely to be a barrier 
transmission artifact for two reasons. First, as we already mentioned above, the maximum 
bias energy is far below the Si02 barrier height, so that finite voltage barrier distortions are 
negligible. Second, this behavior is robust; it has been observed in over 25 junctions on Si:B 



of varying dopant densities, with junction resistances {dV/dl at +1 mV at 1.2 K) ranging 
from 60 kQ to 850 kQ over the same junction area. This order-of-magnitude variation 
in junction resistance is a consequence of uncontrolled variations in barrier thickness and 
purity. The fact that a common high-bias conductance is observed despite such differences 
in barrier properties indicates that the conductance form in Fig. ^ results from the DOS, 
not the barrier. 

Figure § show details of the temperature dependence of the tunneling conductance for the 
86 % sample from 10 K to 0.5 K. Ordinary thermal broadening {i.e. broadening due solely 
to convolution with the df /d{eV) term in Eq. ( p.l|) has not been deconvoluted from these 
data. At 10 K the tunneling conductance is essentially constant within ± 1 mV of V^ = 0. 
As temperature decreases, the opening of a roughly quadratic Coulomb gap is evident. The 
gap has its own temperature dependence (apart from the ordinary thermal broadening), 
becoming both wider and deeper as T decreases. The zero-bias DOS A^(V^ = 0,T) goes 
approximately as T^ at low temperature, turning over and closing {i.e. G{0)/G{+2'mV) ~ 
1) near 7 K. 

III. SCALING ANSATZ FOR THE COULOMB GAP 

Coulomb interactions in a disordered insulator are known to cause a correlation gap in 
N{6). The original ES derivation [^] for the Coulomb gap shape was given for a classical 
disordered system of point-like localized electrons repelling each other via the Coulomb 
interaction U{r) = e'^/nr. This derivation was based on a stability criterion for the ground 
state with respect to transfer of an electron from an occupied state i{ei < 0) at position r^ 
to an empty state j{ej > 0) at position Vj, with: 

e, - £, - f/(r,,) > 0, (3.1) 

where Si and Ej are one electron energies relative to the Fermi level. If states i and j are 
within energy e of the Fermi level, i.e. Ej, \Ei\ < e, then their typical separation in space 
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r{e) = e'^/ne is large He is small. The DOS N{e) ~ d{r{e) ^)/de then has a soft gap around 
e = 0. In 3-D the result is: 

A'W - ~, (3.2) 

In the deep insulator, this quadratic gap extends until the non-interacting value A^o is reached 
(Fig. 1^). The Coulomb gap of Eq. ( |3.2|) leads directly to the ES form for the variable-range 
hopping (VRH) conductivity 

a = aoe-(^°/^)^'^ (3.3) 

where 

To = C-^4' ^3.4) 

and C ^ 2.8 and ^ is the localization length. 

The original ES theory applies to a disordered system far into the insulating side of the 
MIT, However, the vast majority of experimental observations of ES hopping use samples 
doped within 10 % to 50 % of Uc- Therefore, a description is required of how the Coulomb 
gap changes as ric is approached from below. One can formulate a phenomenological scaling 
ansatz which relates the DOS in the insulating state to other critical quantities using the ES 
argument. The two quantities needed to describe both disorder and Coulomb interactions 
near the MIT are ^ and k. In a doped semiconductor, when n ^ Uc from below, ^ and k 
diverge with the decreasing parameter (1 — n/ric) as: 

e(l-n/ne)=a(l-n/n,)-" (3.5) 

k(1 - n/ric) = Ko{i/a)'^'^ = Ko(1 - n/ricY^ (3-6) 

where a = n~^'^ is the average distance between dopants at n = ric, C = ^iv ^ 1); and u 
and 7] are scaling exponents, the same as used by McMillan. |^ (McMillans scaling analysis 
for metals does not calculate explicitly the values of these exponents but does yield the 
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restriction 1 < 77 < 3. Experiments |Q have shown that C ~ 1.) Since k diverges near 
the MIT, the Coulomb interaction becomes weaker and the Coulomb gap becomes steeper 
and narrower. Therefore the VRH conductivity increases due to both the increase of the 
hopping rate at a given distance following Eq. ( |3.5| ), and the increase in the DOS. This is 
reflected in Eq. (|3.4|), where Tq tends to zero as the MIT is approached from below: 

kBToin/ric) ~ A(l - 71/71^^ (3.7) 

where A = e^ / kqu. If r/ ^ 2 and z/ ^ 1, one gets rju ^ 2 which agrees reasonably with DC 
transport measurements in many materials. [33,35| This gives important indirect evidence 



that the DOS is given by Eq. (|3.2| ) with k provided by Eq. (3^). However, DC transport 
is determined entirely by a relatively small range of energies e < ksTo, so that the widest 
possible band of energies which contribute to VRH conductivity is fc^To. Since fcsTo goes 
to zero near ric, conductivity data contain no information about the DOS over a large range 
of energies covering 

ksTo <e <A (3.8) 

In other words, in transport experiments we need only consider energies smaller than /cbTq 
or, equivalently, distances larger than ^. Only at such distances can the hopping probability 
be considered exponentially small and k be considered independent of length scale. 

Tunneling spectroscopy can give direct information about the DOS in the whole range 
e < A, which forces us to account for interactions at distances smaller than ^. This is 
accomplished by introducing spatial dispersion into k at distances r ^ ^ or, in terms of 
wavevector q, at qC, ^ 1: 

K ~ Koir/aY'^-^^ ~ KoiqaY^-'^^ (r < (3.9) 

At r = C, Eq. (|3.9|) matches Eq. (|3.6|) , which is valid for r ^ C^. Eq. ( p.9|) means that at 
r <^: 

U{r) = A{a/ry, (r < (3.10) 
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(One sees from this expression that the index t] is identical to the standard dynamic scahng 
exponent z.) 

Our goal now is to obtain the DOS throughout the energy range ksTo < e < A hj 
repeating the ES argument used to derive Eq. ( p.2| ). For r ^ ^, we can use the potential in 
the form of Eq. ( |3.1(J| ) only if we create wave packets of the size r or smaller from states of 
the size ^. Such a packet costs hD{r)/r'^ of additional "localization" energy, where D{r) is 
the diffusion coefficient at length scale r. We argue below that this energy is smaller than 
the energy gain per particle given by Eq. (p.lO|) and thus allows us to proceed with the ES 



argument. Using Eq. (|3.1| ) and Eq. ( p.lO|) we calculate an average distance r{€) = a{A/eY^^ 
between electron and hole states in the band of energies of width e around ep. Then using 
Eq. (13) we obtain a "critical" DOS: 

N,{e)^a~^A-\e/A)l-\ (3.11) 

At finite ^, Nc{e) matches the quadratic part of Eq. ( |3.2| ) at e = ksTQln). Thus the width 
of the parabolic gap is: 

6 ~ kBToin) ~ A(l - n/n.y' (3.12) 

Right at the transition (where [1 — n/nd = l/C, = 0) Eq. ( ^.11| ) describes the DOS at all 
energies e leg A (which is why we call it critical). Thus at the very transition this scaling 
model has only one unknown index, t]. If t] ^ 2, Eq. ( |3.11|) gives N{e) oc e^^'^. Critical and 



near-critical behavior of the DOS at n = Uc is shown in Fig. |^ by a thick line. This curve 
plays the role of a "backbone" from which N{e) deviates down on the insulating side of the 
transition (n < Uc) and up on the metallic side {n > n^ near e = {i.e. at distances ^ ^). 
These arguments have proceeded from the insulating side. The behavior of the Hartree 
Coulomb gap on the metallic side can also be estimated in the language of the dielectric 
constant, wave packets, and the ES argument. At g,^ ^ 1 (large energies) t here is no 
difference between insulating and metallic phases so that Eq. ( |3.9| ) is also valid for n > ric. 
A qualitative difference between metallic and insulating phases appears only when g^ <^ 1, 
where instead of the scale independent k of Eq. ( p.6|) one gets metallic scaling: |l8|,|36| 
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K{q) ^ Ko{a/0^'-''\qO-' (^^ « 1) (3.13) 

This divergence of K,(q) at small q leads to exponential screening of the interaction at dis- 
tances r > ^: 

f/(r) = (eV«:or)(a/0''-' exp(-r/0, (r « (3-14) 

The role of such screening was studied in Ref. 0. In the first approximation it leads to 
smearing of N{e) at the scale oi 6 = U{^) = ksTQ. Thus in the metallic state N{e) becomes 
finite at e = 0: 

A^(0) = Aa-'^/\-^\n/n^ - l\<^-'i\ {n > nj (3.15) 

The above expressions are all obtained in the Hartree approximation (with excluded self- 
interactions). In this broader sense, all these depletions of the DOS can be called Coulomb 
gaps. However, at very low energies in the metallic state, the Hartree approach fails and 
exchange interactions become dominant, leading to a square-root AA dip in the DOS near 
the Fermi energy. ||^ We can find the limits of the Hartree approximation if we compare the 
energy loss due to creation of small wave packets with the energy gain due to creation of a 
Coulomb gap. The first quantity is hA{r)/r^ and the second is just e. Their ratio: 

Q = hD{r)/r^6 = ha{r)/r^eg{e)e^ = GJr^eg{e) (3.16) 

where a{r) = Gcc'^/h is the conductivity in a scale r, and Gc is the critical dimensionless 
conductance. At the very transition r{e) = a{A/ey^''^, so that Eq. (|3.16| ) with the help of 



Eq. ( p.2|) shows that Q < 1, meaning the Hartree approximation is still reasonable. This 
justifies the backbone DOS Eq. ( p.ll|) at energies e > 6 on both sides of the transition. 



On the other hand, a.t Sc < 6 in the metallic state the Hartree approximation fails, wave 
packets overlap, and an AA-like theory based on exchange interactions becomes valid. As a 
result, a negative square root of comparable amplitude should be added to Eq. ( |3.15| ). 



Now we would like to compare our results with those obtained by McMillan |^. In 
his work the Hartree interaction was neglected and the metallic side of the transition was 

12 



strongly emphasized. He did not present an explicit form of low energy behavior of DOS on 
the insulating side and did not mention that it is described by notion of the Coulomb gap. 
Nevertheless his estimate of the "correlation gap" width A is in agreement with our width 
6 = ksTo of the parabolic Coulomb gap. In the metallic phase and in the critical regime our 
results are identical to those of McMillan |^. (Note that we use same t] and u as McMillan, 
but his A is equivalent to our fcsTo or 6.) 

Refs. P,^ were criticized 0] for using the single-particle DOS N{e) in the dielectric 



constant: 

K = Koil + ^ 3.17 

The standard expression for k of a homogeneous system uses the thermodynamic DOS 
dn/dfi in place of N{e) in Eq. (|3.17|) , where /i is the chemical potential. In an electron 
system with direct Coulomb interactions, dn/dfi is believed to have no gap near the Fermi 
level, in contrast to N{6). So which of the two DOS enters Eq. ( |3.17| ) is the crucial question. 



This question is directly relevant here, since all the results in this paper appear to agree 
with the use of Eq. ( ^.17| ). Indeed, Eqs. (|3.(j|), (^.91), and ( pJ.13D are self-consistently related 



by Eq. ( |gT7[ ) to Eqs. Q, pil| ), and (gJ^ respectively. 

We can repeat the problem in terms of the screening radius Vg by rewriting Eq. (p.l7|) 
as K = kq + l/{rs'^q'^). If r^^ = KQ/ATTN{e)e^ then at small e, N{e) is small, the radius r^ 
is large, the screening is weak so the interaction is strong, and finally N{e) has a Coulomb 
gap. For example, deep on the insulating side, using Eq. ( |3.2|) we get that at energy e the 
radius r^ = e^ / kqE, exactly equal to the average distance between states within e of the 
Fermi level. This implies that screening self-consistently does not destroy the Coulomb gap. 



The other option argued for in Ref. |T0| starts from the DOS dn/dfi. If dn/dfi is large, 
^s^ = Ko/(4vre^(ira/(iyu) is small and we get a short range interaction which does not lead to 
the Coulomb gap. This result is inconsistent with the ES argument, computer simulations, 
and numerous observations of the ES law for VRH conductivity. Why, then, does the 
conventional definition of the screening radius fail? 
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An answer to this question for the deeply insulating phase was suggested in Ref. [^ . 
While the correct expression for n contains dn/dfi, in a strongly inhomogeneous system with 
localized states one needs a local k, and dn/dfi to describe the interaction of two particular 
states at the distance r. The local {dn/dfi)~^ fluctuates strongly in absolute value and has 
random sign. The inverse amplitude of these fluctuations equals N{e). Thus, the local 
screening is determined by a small random sign DOS which in absolute value is of order 
N{e). The random sign of the interaction does not change any estimates based on the ES 
argument, and this is why the ES argument works. On the other hand, {dn/dfi)~^ averaged 
over realizations of a random potential is very small, so that the thermodynamic DOS 
dn/dfi is large and energy independent. The absence of self averaging in strongly disordered 
systems is the main reason that the average dn/dfi does not determine the strength of local 



interactions. Ref. |^ discussed a classical case of a deep insulator. We assume that in the 
critical range of the MIT quantum effects can be taken into account by renormalization of 
^ and K, and by introduction of wave packets. After that we have a problem similar to the 
classical one, but with renormalized interactions. This interaction for a given realization of 
impurities is not screened. To find N{e) one needs the local interaction, and the average 
dn/dfi is irrelevant. 

Let us remind the computer simulations lead to these conclusions [^. In that work the 
authors considered randomly distributed pointlike donors and acceptors, all acceptors being 
negatively charged and some donors being occupied by electrons (and hence neutral) and 
others empty (and so positive). The ground state of such system is a classical realization of 
a Coulomb gap [0,0. Screening a probe positive charge at T = was studied numerically. 
It was shown that for a given realization of impurities, the screened potential of the probe 
charge is almost random in sign and its absolute value decays as 1/r. To understand the 
origin of the random sign potential for a single realization of potential, recall that in the 
ground state electrons are correlated in space or, in other words, positive donors alternate 
with neutral ones. The probe positive charge can attract a new electron to the nearest empty 
donor. The positive donor left behind attracts a new electron to another of its neighbors, and 
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so on. As a result a chain of new alternating positive and negative charges can appear. Such 
a chain has random direction and length and no symmetry. It is clear that the sign of the 
potential, for example at the end of the chain, is random. Huge cancellations appear when 
we average over many different realizations of impurities or over different configurations 
of chains. As a result spherical symmetry is restored and an exponentially small positive 
potential survives. This picture has been generalized for a description of the Rudderman- 



Kittel exchange interaction in a disordered metal. [^ 



An alternative explanation was recently suggested by Si and Varma, |^ who argued 



that, contrary to previous assertions, direct Coulomb interactions do lead to a gap in the 
average dn/dfi when the density is low enough in a 2D system. They suggest that the 
incorporation of this result into the charge diffusion models of Refs. |[T^-[T^ could generate 
a realistic description of the transition from diffusive metal to Coulomb glass insulator. 

Finally, Fig. ^ presents numerical calculations of DOS spectra of the impurity band model 
obtained by computer minimizations of the total energy of the system in the framework of 
our scaling ansatz model. To describe the MIT at n < ric, n = Uc, and n > Uc, we use 
Coulomb potentials modified by Eqs. (PTB|), (^.91), and ( |3.13| ) respectively. For concreteness. 



we assume i] = 2 and f/(r) = (R/ p^ for n = ric, U{r) = (p^ + ^'^y/'^R/ p^ for n < n^; 
f/(r) = i?^/[p^(expp/^ — 1)] for n > ric Here R = n~^^^ is the average distance between 
donors at the transition, U is measured in units of c'^/kqR, and p^ = r"^ + (i?/2)^. The 
(i?/2)^ term in the last expression is significant only at small r and is introduced as a large 
energy cut-off (or alternatively as a small wavelength cut-off at distances a ~ R/2 at Uc). 
The algorithm of the simulation is described in Ref. [0]. The energy reference is adjusted 
for each realization of impurity coordinates so that the chemical potential is zero, and the 
DOS is averaged over 10000 realizations. In Fig. ^ the insulating and metallic DOS are 
computed for states equidistant from the transition, using ^ = 2R. We see that results of 
calculation are in a reasonable agreement with the schematic scaling plot of Fig. ^ and the 
experimental data of the previous section. 
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IV. ANALYSIS OF THE TUNNELING DATA 

There are clearly remarkable similarities between the experimental data of Fig. ^ and the 
scaling ansatz results in Fig. ^. In fact, most major features of the tunneling measurements 
can be accounted for by the simple scaling analysis of the preceding section. It is obvious from 
Fig. ^ that the high-bias conductance common to both localized insulator and disordered 
metal plays the role of the high energy backbone DOS obtained from the scaling ansatz 



by combining the dielectric constant of Eq. (|3.9|) with the Coulomb gap of Eq. (p.2|). The 



distinction between metal and insulator becomes clear only when the DOS departs from 
the backbone at low energies. From the data, this departure occurs near an energy scale of 
about 0.5 meV. 

While the scaling ansatz cannot calculate the value of the exponent describing the power 
law backbone DOS, the experimental data does yield this number and the value of the 
associated scaling exponent. According to Eq. ( p.lO|) , the high-energy backbone DOS near 
the critical density for both metals and insulators depends on a single exponent rj and should 
go as Nc{e) oc e^~ . From the data of Fig. ^ in the bias range 1 mV < V < 50 mV we 
obtain G{V) oc V"^ with values for the exponent m = 0.45 ± 0.2 covering all the samples 
measured. This gives rj = 3/(1 + m) = 2.1, which is certainly within the theoretically 
required bounds 1 < ?7 < 3 and in fact is essentially the same as the value 77 = 2 reported by 



Hertel, etal. ||TI| in barely metallic NbSi alloys. Since this value for 77 was obtained from the 
high-energy DOS, where metals and insulators share a common DOS character, and from 
many different samples both below and above ric, we believe it is a reliable value. We do 
wish to emphasize that, in our interpretation, the approximately square-root behavior of the 
DOS at these high energies is unrelated to the AA exchange correction, which is important 
only at energies e < 6 ^ 0.4 meV in the metallic state. 

The exponent z/ can also be examined, though less definitively than 77. On the metallic 
side of the MIT this can be done using Eq. (|3.15| ). The measured ratio of zero-bias conduc- 
tances for the two metallic samples G'iio(0)/G'io5(0) = 1.7. Relating this conductance ratio 
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to the ratio of the DOS and using Eq. ( p.l5|) , we obtain z/(3 — 77) ^ 0.77. Taking rj ^ 2 from 
the preceding analysis, this gives z/ ^ 0.77. Ahernatively, z/ can be estimated independently 
using data from the insulating samples. Fig. |^ shows the measured Coulomb gap widths 
for four insulating samples plotted against normalized dopant density. The data fit a linear 
function very well with nearly zero intercept, i.e. 6{n) oc (1 — n/ric)'^ where 7 ^ 1.0. From 
Eq. (|3.12|) , we have 7 = rju, so that if 77 = 2, then the gap width data imply that u = 0.5. 
Fig. also shows that, in agreement with Eq. (|3.12|) , the values of Tq measured from DC 



conductivity in the more insulating samples, where the ES hopping law is followed most 
rigorously, are of the same order as the corresponding values of 6. 

These values for u obtained from the tunneling data can be compared to the independent 
data for K,{n) shown in Fig. |l]. A simple power-law fit to the data in that figure gives 
( ~ 0.71 = z/(?7 — 1). Again, if 77 is taken to be 2, then u ^ 0.71. This value is comparable 
to the values obtained from the DOS data analysis. Therefore, we conclude that the DOS 
data support values for the scaling exponents of r^ ^ 2 and z/ somewhere between 0.5 and 
0.8, at least within the range of n/uc covered by the samples used. 

The exponent z/ is commonly measured by examining the very low temperature DC con- 
ductivity a(T — * 0) as a function of doping and fitting the data to the form cr(0) ~ {n — ricY 
for n slightly above Uc. For nominally uncompensated doped silicon (Si:B, 1^^ Si:P, ^ 



and Si:As [^), many such transport experiments have reported values for u between 0.5 
and 1, generally closer to 0.5 than to the u ^ 1 observed in most other disordered conduc- 
tors and expected theoretically. In Si:B, for example, the exhaustive very low temperature 
conductivity measurements reported in Ref. [^ indicate that z/ = 0.65. Thus our analysis 



of the tunneling data yield a value for u at least consistent with that obtained from trans- 
port measurements. However, the interpretation of the transport data leading to z/ ^ 0.5 



in Refs. p7| , ^ , ^ has recently been questioned. P^P3[] There is some evidence that if only 
samples within ~ 1 % of nc are used, then transport measurements yield z/ ^ 1 instead. In 
fact, stress tuning of a Si:B sample ^^ across the MIT has been reported to give u as high 



as 1.6, although earlier stress experiments [^ on Si:P reported u = 0.5. 
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An estimate for the value of ^ can also be extracted from the data. Using Eq. ( p.4|) 
combined with the measured values of /t and Tq in Figs. |I| and |^ respectively, we obtain ^ ^ 
25 nni, 30 nni, and 40 nm for the 81 %, 86 %, and 90 % samples, respectively. These values 
are in reasonable agreement with Eq. ( |3.5| ) using v ~ 0.7 and length a = 7.5 nm, which is 
close to the average distance between impurities at Uc, 6.3 nm. 

It is necessary to discuss whether the scaling ansatz of Section |T| should be really ap- 
plicable to the experimental situation of Section ||. First there is a problem with theweak 
compensation of the Si:B samples used. Indeed the theoretical discussion starts from the 
picture of a classical Coulomb gap. In the impurity band of a lightly doped semiconductor, 
the Coulomb gap is most pronounced for degrees of compensation (ratio of minority to ma- 
jority impurity concentration) comparable to 1/2. (In this case a large number of impurities 
are charged, but the peak of the DOS is still not too broad. 0) Then with increasing n 
the Coulomb gap evolves according to the predictions of Sec. 3. In uncompensated lightly 
doped semiconductors, all majority impurities are occupied in the ground state and there- 
fore are neutral. There is a large Mott-Hubbard gap which separates occupied states from 
empty states of a second electron on an impurity. Disorder is very weak because there are 
no random charges. So the Mott-Hubbard gap at the Fermi level of a lightly doped un- 
compensated semiconductor is a real hard gap as opposed to the soft Coulomb gap. Why 
then can we then talk about a Coulomb gap near the MIT? Bhatt and Rice |^ pointed 
out that at large impurity concentrations n ~ rzc a new phenomenon which they called 
self- compensation can take place. Because of large overlap and strong positional disorder 
of their wave functions, some clusters of impurities can have large affinity to electrons and 
can ionize other clusters. As a result charges will appear and a random potential will close 
the Mott-Hubbard gap which is already narrowed at these concentrations. Then our theory, 
developed for a compensated semiconductor can work even for a nominally uncompensated 
one. Similar tunneling experiments on compensated samples would be very important to 
find out whether they provide a similar DOS and whether it agrees with our theory. 

A second problem is related to the role of the Al electrode in the screening of the long- 



range Coulomb interactions. One can think that tunnehng electrons typically penetrate 
a distance ^ into a semiconductor. The characteristic length r{e) of interactions which 
determines the DOS at e ^ 6 (square root region) is smaller than ^. This means that the 
large square root region is robust with respect to screening by a metal gate. However, the 
parabolic gap results from the Hartree interaction between charges separated by distances 
r > ^. Therefore, the Al electrode could produce substantial image charge screening of this 
interaction. Such screening was studied numerically [Q for the case of a classical Coulomb 



glass and a metal electrode placed directly on its surface. These computations show the 
bulk parabolic Coulomb gap near the surface can be closer to a linear one, rather than 
quadratic one. In our experiment, electrode screening may be reduced by two factors. First, 
the Al electrode is backed off the semiconductor surface by roughly 1.5 to 3 nm of Si02, as 
described in Section 0. Second, the Schottky barrier between the Al and Si can be as wide 
as 20 nm. The cumulative effect of these two barriers is to increase the distance between 
the charges and image charges to distances ~ ^. However, even if the gate is backed off, 
the parabolic gap should still be expected to fail at small enough energies. Surprisingly, 
the parabolic gap is more robust than expected. There is strong agreement between the 
power law for the Coulomb gap as measured by tunneling and by bulk DC transport, which 
is not affected by electrode screening because the sample is macroscopically thick (~ 250 
to 300 yum) so that most of current flows far from the Al contacts. If the Coulomb gap is 
described by N{e) ~ e^, then the ES model gives a hopping exponent s = {p + l)/{p + 4) 
for the conductivity a = (Tq exp{ — (Tq/T)'^}. For the 86 % sample at 0.5 K, a power-law 
fit to the tunneling data gives p = 2.2, so the predicted ES value for s = 0.52. Transport 
measurements on the same sample in the same temperature range give a value s = 0.51±0.02. 
The close consistency o the tunneling and transport results indicates that, to a reasonable 
accuracy, the tunneling DOS reflects the bulk characteristics. 

Finally, we should mention that the definition of DOS discussed in section 2 does not 
strictly coincide with the definition of the tunneling DOS measured experimentally. Sec- 
tion Jl| dealt with the bulk DOS of charged excitations (electronic polarons). For a deep 
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insulator these excitations were introduced in Ref. [^. Electronic polarons are responsible 
for ES VRH. On the other hand tunneling experiments measure the one-electron DOS. For a 
deep insulator, a substantial difference between these two DOS in the limit of small energies 
was predicted in Ref. [0 , but was not reliably observed in numerical simulations. We cannot 
resolve here this complicated problem. It may happen that this difference can explain the 
apparent absence of the effect of the gate on the parabolic gap. We would like only to 
emphasize here that this uncertainty has to do only with small energies e -C 5 inside the 
parabolic Coulomb gap. 

V. SUMMARY 

Electron tunneling measurements of the single-particle DOS have beenmade on Si:B 
crystals ranging from 81 % to 110 % of the critical density of the MIT. At low energies {e < 
0.5 meV), non-metallic samples show an approximately quadratic soft Coulomb gap across 
the Fermi level, while metallic samples show a square-root cusp. At higher energies, (up to 50 
meV), both insulating and metallic samples show a common DOS behavior, N{6) ~ e'"^. with 
m slightly less than 1/2. These features of the data can be understood within the framework 
of a scaling ansatz of the approach to the MIT from the insulating side. By extending the 
semi-classical Efros-Shklovskii derivation of the Hartree Coulomb gap to include a spatial 
dispersion of the dielectric constant at small length scales, it was shown that at short length 
scales (r <^ ^), or high energies, spatial dispersion of k gives rise to a common DOS that 
increases as a power law in both metals and insulators. This is distinct from the Altshuler- 
Aronov low-energy square-root cusp in the DOS arising from exchange correlations in weakly 
disordered metals. Only at long length scales, or, equivalently, low energies, does a quadratic 
Coulomb gap become apparent in insulators, while metals take on a non-zero DOS across 
the Fermi level. Finally, the tunneling results and the scaling ansatz suggest that a semi- 
classical approach may be used to treat Coulomb interaction effects on the approach to the 
MIT from the insulating state. 
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FIGURES 
FIG. 1. Static dielectric constants (in units of the Si dielectric constant kq = 11-7) as a function 

of the normalized dopant density for a set of insulating Si:B crystals. These data were measured at 

a temperature T = 0.1 K. The solid curve is a fit to the functional form K{n/nc) -^ (1 — n/nc)~, 

which gives C ~ 0.71. 

FIG. 2. Low-bias tunneling conductance-voltage spectra for Si:B samples at several different 
values of n/ric- All data sets are normalized to the conductance at +2 mV. Data for samples with 
n/ric > 90 % were taken at T = 0.1 K. Data for the 86 % and 81 % samples were taken at 0.5 K 
and 1.2 K, respectively. 

FIG. 3. Extended bias range tunneling conductance spectra for four Si:B samples, two metallic 
(110 % and 105 %) and two insulating (95 % and 90 %), all taken at temperature T = 0.1 K. The 
data are plotted on a log-log scale and are normalized to the conductance at +50 mV. 

FIG. 4. Low-bias tunneling conductance-voltage spectra for the 86 % Si:B sample at several 
different temperatures. The spectrum at each temperature is normalized to the conductance at +2 
mV. Ordinary thermal broadening has not been deconvoluted. 

FIG. 5. Schematic plot of the DOS as a function of energy for different n: 1.) n < nc/2, 2.) 
n < Uc and ric — n <^ ric, 3.) n = ric, 4.) n > ric and n — ric <ti ric, and 5.) n > nc/2. 

FIG. 6. DOS as a function of energy as obtained by numerical simulation for the Coulomb 
potential (dashed line) and for three model potentials corresponding to the critical range of the 
MIT (full lines) from bottom to top: n < nc{^/R = 2), n = Uc, and n > nd^/R = 2). 

FIG. 7. Full width 6 of the parabolic Coulomb gap in four insulating samples plotted against 
(1 — n/ric). The dashed line is a linear fit. 



26 



20 
18 
16 

v° 14 

y 

12 
10 



1 

□ 


1 1 
Si:B 


- 


T = 0.1K 


_ 


T\ 1 (l-n/n)-"-^' 


1 


'■+ , 

1 1 i 



0.1 
1 - (n/n ) 



0.2 



l-O^I^ 1 1 1 1 1 1 1 _|i 
^_ n/n = 110% m.^ 


1.0 


1 1 1 1 1 1 1 «< 

.to 
•to - n/n = 90% »•• 
- •to. •• 


0.8 


•• «• 


0.8 




• T = 0.1 K •• 




•to. T = 0.1K^to* 


0.6 


\ .•• 


0.6 


0.4 


V 


0.4 


to to 

to to 
to to 


0.2 




0.2 


to to 





1 1 1 1 1 1 1 





^^^ 


1 




k 1 




< 
0.8 


. 1 1 1 1 1 1 \ d^ 
•• _ n/n = 105% • • 

•a c # _ 

\ T = 0.1 K •• 


' i'J 1 1 1 1 1 1 1 ^W" 

••• 

0.8 -••to n/n^ = 86% ,• 


6 


-.. 


0.6 


•» T = 0.5 K • 

•to to 




\ / 




1 1 

to 

1 1 


0.4 


V 


0.4 


0.2 


• 


0.2 


0.0 


1 1 1 1 1 1 1 


0.0 


^r 


1.0( 
0.8 


(■^ 1 1 1 1 1 1 1 ^^ 
%y^ n/n^ = 95% ^f^ 
^W T = 0.1K ^F 


. 1.0 
0.8 


1 1 1 1 1 1 '^■l'^ 
^Nto. n/n =81% »• 

^m • 

• T=1.2K • 


0.6 


\ / 


0.6 


to to 

to to 
to to 


0.4 


\ / 


0.4 


to to " 




• • 




to to 


0.2 


\2 


0.2 


to to 




• • 




"•toatoto* 





\ tm \ 





1 1 1 1 1 1 1 


-2-10 1 1 


i -2 -1 1 2 




V (mV) 




V (mV) 



> 



+ 



> 



1.0 
0.8 
0.6 
0.4 
0.2 



0.0 



-2.0 -1.5 




o 


10 K 


♦ 


4.2 K 


n 


1.0 K 


• 


0.5 K 



-1.0 -0.5 0.0 0.5 
V (mV) 



0.2 



^0.1 



0.0 

-0.4 



\ 1 




1 






1 y 


\ \ \ \ 
\ \ \ \ 










/ y\^ 


\ \^ ^\ \ 












\ \ \ \ 
\ \ \ \ 
\ \ \ \ 
\ \ \ \ 


\ 






/ 


yy,/' 




\ 






/ 


/ / / 


\ \ \ 












\ \ \ 












\ \ \ 






















/ y / 


\ \ X 








/ 


/ / 


■^ \ ^ 








/ 


/ / 


X \ 

\ 

X \ 


\ 






/ 


/ 

/ y 
y 
y 


^^ 












\ 




X. \ 






y^ 


1 


^x^ 


.^Vi 


^ — -^ 


y 


1 



-0.2 



0.0 

e 



0.2 



0.4 



2.5 
2.0 



> 



u 
c 
W 



1.0 



0.5 - 



0.0 
0.00 





1 
Si:B 




1 

1 


- 






.-^- 










.^^ ^ 


-A -5 












I'- 


- 


• k T 

B 




y 


\ 






I-" 




" 


y 


\ 


» 




r 
1 


. 


1 



0.05 0.10 0.15 

1 - (n/n ) 



0.20 



